What drives mesoscale atmospheric turbulence? 



H. Xia, 1 H. Pimzmann, 1 G. Falkovich, 2 and M.G. Shats 1 ^ 

1 Research School of Physical Sciences and Engineering, 
Australian National University, Canberra ACT 0200, Australia 
2 Physics of Complex Systems, Weizmann Institute of Science, Rehovot 76100, Israel 

(Dated: May 4, 2008) 

Measurements of atmospheric winds in the mesoscale range (10-500 km) reveal remarkably uni- 
versal spectra with the k~ 5 ^ 3 power law. Despite initial expectations of the inverse energy cascade, 
as in two-dimensional (2D) turbulence, measurements of the third velocity moment in atmosphere, 
suggested a direct energy cascade. Here we propose a possible solution to this controversy by ac- 
counting for the presence of a large-scale coherent flow, or a spectral condensate. We present new 
experimental laboratory data and show that the presence of a large-scale shear flow modifies the 
third-order velocity moment in spectrally condensed 2D turbulence, making it, in some conditions, 
similar to that observed in the atmosphere. 
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Atmospheric motions are powered by gradients of so- 
lar heating. Vertical gradients cause thermal convec- 
tion on the scale of the troposphere depth (less than 10 
km). Horizontal gradients excite motions on a planetary 
(10000 km) and smaller scales. Both inputs are redis- 
tributed over wide spectral intervals by nonlinear inter- 
actions [U El E]. The spectra of kinetic energy of at- 
mospheric winds have been analyzed during the Global 
Atmospheric Sampling Program [3J. These wavenumbcr 
spectra measured in the upper troposphere and in the 
lower stratosphere have shown two power laws: E(k) oc 
f or t ne s Ca les between 10 and 500 km, and a steeper 
spectrum with E(k) oc fc~ 3 in the range of scales (500- 
3000) km (similar to the spectra in Figs. 2,3). Ac- 
cording to the Kraichnan theory of 2D turbulence [JJ, 
E(k) = Ce 2 / 3 fc -5 / 3 corresponds to an inverse energy cas- 
cade and E(k) — C q r) 2 ^ 3 k~ 3 to a direct vorticity cas- 
cade. Here e and rj are the dissipation rates of energy 
and enstrophy respectively. This prompted a two-source 
picture of atmospheric turbulence with a planetary-scale 
source of vorticity and depth-scale source of energy [5]. 
The large-scale part of the spectrum can be due to a di- 
rect vorticity cascade [5] or it can result from an inverse 
cascade of inertio-gravity waves [B]. It can also result 
from the energy pile-up at the system scale in the pro- 
cess of spectral condensation [7] (a turbulent counterpart 
of Bose condensation) , or it can be due to a combination 
of the above. What about the mesoscale 5/3-spectrum? 
Is it an energy cascade and what is the flux direction? 

In homogeneous turbulence, spectral energy flux is ex- 
pressed via the third-order moment of the velocity [8j |9] 
S3 = ((SV L ) 3 ) + (SV L (SV T ) 2 ) . Here 5V L and SV T denote 
the difference of velocities at two points separated by 
distance r. Angular brackets denote ensemble averaging 
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over realizations, and the subscripts denote the longitu- 
dinal (L) and transverse (T) velocity components rela- 
tive to r. Positive S3 corresponds to the inverse energy 
cascade from small to large scales. Measurements of the 
third moment of the velocity difference in the atmosphere 
gave a negative value in the interval 10-100 km, which 
was interpreted as the signature of the forward cascade 
|10j . The negativity of the third moment spawned hy- 
potheses about a direct energy cascade in 2D or stratified 
turbulence [TTJ [T5] . Here we demonstrate experimentally 
that a negative small-scale S3 in a system with an inverse 
cascade can be caused by a large-scale shear flow. 

Let us first consider how small- and large-scale parts 
of the velocity difference (respectively Sv and SV) con- 
tribute to the second and third velocity moments. Com- 
paring ((SV) 2 ) = s 2 r 2 with ((5v) 2 ) = C(er) 2 / 3 we see 
that the small-scale (turbulent) part dominates at the 
scales smaller than l t = C 3 ' 4 s — ^e 1 / 2 . Here s = V/L s is 
a large-scale velocity gradient and L s is the velocity shear 
scale length which depends on the system size and on the 
topology of the large-scale flow. For the third moment, 
we compare ((Sv) 3 ) = er with the cross-correlation term 
(SV(Sv) 2 ) = srC(er) 2 / 3 and observe that the influence of 
SV extends to a much smaller scale Z* = C~ 3 / 2 s~ 3 l 2 e 1 / 2 , 
because the dimensionless constant C is typically larger 
than unity, as discussed below. 

The above estimates are true for a large-scale part pro- 
duced by any source. In particular, when it is produced 
as a condensate by an inverse cascade [13j [T4J [15l [16j 
[17] 118] (as in the experiments described below) one es- 
timates s as follows. Let the linear damping rate a be 
smaller than the inverse turn-over time C 1 / 2 e 1 / 3 L~ 3 / 2 
for the vortices comparable to the system size L. Then 
the flow coherent over the system size (the condensate) 
appears [THJ El HS1 dl H3 HB] with the velocity esti- 
mated from the energy balance, aV 2 s» 2e, which gives 
s = VIL. = Lz x -y/2du and 
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Note that this is not the condition that the turnover time 
at It is a" 1 , as in [13]; incidentally Eq. [l] gives a correct 
estimate (k t = 1) for their conditions. The spectrum 
E(k) oc k~ 3 at k < k t is due to the condensate [T3J [T5J 
US], while = Ce 2 / 3 /;" 5 / 3 is expected at fc t < k < k f . 

Here we report the experiment in which the strength 
and the spectral extent of the condensate are varied by 
changing either a or L. The experimental setup shown in 
Figure [T] is similar to those described in [T7J [TSJ but 
has a substantially larger number of forcing vortices (up 
to 900), higher spatial resolution and larger scale separa- 
tion (L/lf « 30). A turbulent flow is generated electro- 
magnetically in stratified thin fluid layers whose thick- 
nesses are varied to achieve different damping rates a. A 
heavier non-conducting fluid (Fluorinert, specific gravity 
SG =1.8) is placed at the bottom. A lighter conduct- 
ing fluid, NaCl water solution (SG =1.03), is placed on 
top. Non-uniform magnetic field is produced by a square 
matrix of 30 x 30 permanent magnets (10 mm apart). 
The electric current flowing through the top (conducting) 
layer produces 900 J x £?-driven vortices which generate 
turbulence. Square boundaries with L = (0.09 — 0.24) 
m are used. To visualize the flow, imaging particles 
(polyamid, 50/xm, SG= 1.03) are suspended in the top 
layer and are illuminated by a 1 mm laser sheet parallel 
to the fluid surface. Laser light scattered by the particles 
is filmed from above using a 12 Mpixel camera. Green 
and blue lasers (A = 532 nm and A = 473 nm) are pulsed 
for 20 ms consecutively with a delay of (20-150) ms. In 
each camera frame, two laser pulses produce a pair of 
images (green and blue) for each particle. The frame 
images are then split into a pair of images according to 
the colour. The velocity fields are obtained from these 
pairs of images using the cross-correlation particle image 
vclocimetry technique. The velocity fields are measured 
every 0.33 s (at the camera shooting rate). For a bet- 
ter time resolution a video camera (2 megapixel) with a 
single laser is used. The damping rate (in the range of 
a = (0.05 — 0.5) s _1 ) is estimated from the decay of the 
total kinetic energy E after switching off the forcing at 
t = to: Et = ^ a e- a C*-*o). 

Figure [2] shows the energy spectrum measured in the 
large box L = 0.235 m at an intermediate damping of 
a = 0.16 s _1 . A forcing scale corresponds to kt « 400 
m -1 . At k > kf , E(k) cx fc~ 3 , while at k t < k < kf, 
E(k) cx fc~ 5 / 3 . At k < k t ~ 80 m _1 , in the condensate 
range, the spectrum is steeper and close to fc~ 3 . Due 
to the presence of the condensate, the spectrum has k~ 3 
and fc -5 / 3 ranges for the large and intermediate scales re- 
spectively, similarly to the Nastrom-Gage spectrum [4]. 
Spectra for different L and a are shown in Fig. [3] At 
fixed a = 0.15 s _1 , the knee of the spectrum shifts from 
k t ~ 80 m" 1 for L = 0.235 m to k t ~ 135 m" 1 for 




FIG. 1: Schematic of the experimental setup 
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FIG. 2: Kinetic energy spectrum measured for the largest box 
L — 0.235 m and intermediate damping a — 0.16s _1 . The 
guide lines show the power laws for different spectral ranges: 
the k~ 3 vorticity cascade, the k~ 5 ^ 3 energy cascade and the 
fc~ 3 condensate range. 

L = 0.15 m, Fig. [3ja). When L is constant, linear 
damping affects kt as shown in Fig. [3jb) . Going from 
a = 0.15s -1 to a = 0.06s -1 , changes k t from 80 to -1 
to kt ~ 130 m _1 . These observations are in good quali- 
tative agreement with Eq. [l] By further reducing L, we 
can achieve a regime when k t w kt, and the fc -5 / 3 range 
disappears, such that the entire spectrum is E(k) oc k~ 3 , 
both above and below kf, as in [5D]. The closer kt is to 
kf, the more symmetric the coherent flow becomes (e.g. 
circular [20]). Therefore, we can control the shape of 
the spectrum and the relative strength of the condensate 
with respect to turbulence. 

We now analyse two regimes whose spectra are shown 
in Fig. |3|a,c). A weaker condensate of Fig. |4|a) was 
generated at a = 0.3 s -1 and L — 0.235 m, while a 
stronger condensate (Fig.[4jc)) was obtained at a = 0.15 



3 



10 



10 



10 



10 



\E(k) 






1=0.15 m 


1=0.235 m V 




k 






FIG. 3: Kinetic energy spectra (a) for different box sizes at 
a — 0.15 s _1 , and (b) for different linear damping rates at 
L = 0.235 m. 

s -1 and L = 0.15 m. The third-order velocity moments 
differ markedly for these two cases (Figs.[4]jb,d)). For the 
weaker condensate case (Fig. |4|b)), S 3 is a linear func- 
tion of the separation distance r in the range of scales 
corresponding to the k~ 5 / 3 spectral range. Both longitu- 
dinal S^l and transverse S$r moments are positive in the 
entire range of scales, in agreement with expectations for 
the inverse energy cascade. The slope of S3 (r) changes at 
about r = 0.04 m. This scale corresponds to the knee in 
the energy spectrum at kt ~ 80 m _1 . The Kolmogorov 
constant can be determined as C = E(k)e~ 2 ^ 3 k 5 ^ 3 , where 
e = (S 3L + S 3T )/(2r). For the data of Fig. Qb) at 
r < n/kt ~ 0.04 m, we have C = 5.6, which is close 
to the value C = (5.8 — 7) previously obtained in numer- 
ical simulations of 2D turbulence (see 21J and references 
therein). At larger separations, r > n/kt, the function 
S 3 (r) grows faster and deviates from linear. 

For the stronger condensate, the spectrum scales as 
E(k) oc fc~ 5 / 3 in the range kf > k > k t ~ 125 m _1 . 
In this case, S 3 (r) changes sign at r = n/k t (Fig. Qd)). 
Such Ss(r) dependence resembles the third-order struc- 
ture function measured in the lower stratosphere [TO] , 
Note that in our case all the driving comes from small 
scales and there is no direct cascade at all, yet S3 is 
strongly modified compared with the weak condensate 
case. 

The spectral condensation can be viewed as the gen- 
eration of the mean flow, which can be revealed by a 
temporal averaging of the instantaneous velocity fields: 
V(x,y) = ^/NJ2 n =i V(x,y,t n ). The power spectrum of 
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FIG. 4: Energy spectra E^(k) and the 3rd-order velocity 
moments Ss(r) of turbulence with a weaker (a,b) and stronger 
(c,d) condensates. After subtracting mean flow in the stronger 
condensate case: kinetic energy spectrum (e) and 3rd-order 
velocity moment (f). 

the mean flow [T5J [TB] is close to E(k) oc k~ 3 . The ve- 
locity field contains both the mean component and tur- 
bulent velocity fluctuations: SV — SV + SV. From the 
data of Fig. Qc) (stronger condensate) , we estimate that 
((SV) 2 } differs from ((SV) 2 ) by about 20-30%. However 
((SV) 3 ) and ((^V") 3 ) differ by orders of magnitude and 
even the sign can be different. It should be noted that 
signs, values and functional dependencies S 3 (r) vary a 
lot for different topologies of the condensate flows (e.g., 
dipolc, or monopolc) and also depend on the mean shear 
in such a flow. We observe negative and positive S^r) 
in the entire range of scales, as well as the sign reversal, 
as in Fig. Qd). 

To recover the statistical moments of the turbulent ve- 
locity fluctuations we take N = 350 instantaneous ve- 
locity fields, subtract their mean flow and then compute 
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the Fourier spectrum and the structure functions. The 
result for the stronger condensate of Fig. |4jc,d) is shown 
in Figs. |4je,f). The subtraction of the mean restores the 
£-5/3 ran g e j n f ac t ; the E(k) scatter in Fig.^e) is less 
than in the total spectrum of Fig. |4jc) . The subtraction 
has even more dramatic effect on S3. As seen in Fig.Qf), 
S3 is a linear function of r in the "turbulence" range. The 
spectral energy flux is deduced as e = S3/7*. The value 
of the Kolmogorov constant C = i?(fc)e~ 2 / 3 fc 5 / 3 « 7.7 
appears to be slightly higher than in the weak conden- 
sate case, but is still close to the values obtained in nu- 
merical simulations. At the border between condensate 
and turbulence, r t w n/kt, the S^r) dependence changes 
radically, suggesting that the condensate affects the in- 
verse energy cascade at large scales. The recovery of 
the linear positive Ss(r) has also been observed at even 
stronger condensates (e.g. circular monopole flows). The 
Kolmogorov constant seems to become smaller at higher 
flow velocities and stronger shear, C = (1 — 4). 

It is important to note that similarity of our spectra to 
those of [21 H] does not necessarily mean that k~ 3 spec- 
trum at large scales in the Earth atmosphere is also fed 
by the inverse cascade. To establish whether this is the 
case, one needs to analyze the atmospheric data in the 
way described here: subtract the coherent flow, recalcu- 
late the second and the third moment of fluctuations and 
use Eq. [l] It is likely that the baroclinic (large-scale) in- 
stabilities play a role in forcing the large-scale flows. To 
model an extra large-scale forcing in our experiments we 
added a large magnet on top of the small-scale forcing 
(as described in [IS]) and found that the modifications 
in 53 are similar to those when the large-scale flow is 
formed via spectral condensation. The mean subtrac- 
tion recovers the energy flux from small to large scales 
in both cases. Similarly, the mesoscale turbulence in the 
Earth atmosphere should be affected by the large-scale 
flow regardless of its origin. 

Recent numerical simulations (see [12] and references 
therein) show that stratification may enforce a 3D dy- 
namics and the forward energy cascade. On the other 
hand, recent experimental studies of decaying turbulence 
suggest a strong role of rotation in establishing a quasi- 
2D regime in which geostrophic dynamics is dominant 
(regime of low Froude and Rossby numbers) [25]. More 
experiments in constantly forced turbulence are needed 
to better understand the competing effects of rotation 
and stratification along with the complex interplay be- 
tween turbulence and waves, resonant wave- wave inter- 
actions, etc. An ultimate answer to the question asked in 
the title of this Letter can only be resolved in the atmo- 
spheric measurements. What we have shown here is the 
need to separate mean flows and fluctuations to recover 
the energy flux. 
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